Astronomy & Astrophysics manuscript no. reggiani'rn' final 
April 24, 2013 



©ESO2013 



Universality of the companion mass-ratio distribution 

(Research Note) 

M. Reggianff and M. R. MeyefU 



o 

(N 

< 

m 

(N 



Institute of Astronomy, ETH Zurich, CH-8093 Zurich, Switzerland 

Preprint online version: April 24, 2013 

ABSTRACT 

Aims. We present new results regarding the companion mass-ratio distribution (CMRD) of stars, as a follow-up of our previous work. 
Methods. We used a maximum-likelihood-estimation method to re-derive the field CMRD power law avoiding dependence on the 
arbitrary binning. We also considered two new surveys of multiples in the field for solar- type stars and M dwarfs to test the universality 
of the CMRD. 

Results. We found no significant differences in the CMRD for M dwarfs and solar-type stars compared with previous results over the 
common mass ratio and separation range. The new best-fit power law of the CMRD in the field, combining two previous sets of data, 
is dN/dq oc (f, with/? = 0.25 + 0.29. 
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A large portion of stars, both in the field (Raghavan et al. 2010; 
Jans on et alj|2012l) and in star- forming regions (Patience et alj 
2002), are formed in multiple systems. Therefore understanding 
multiple star fo rmation is necessary to investigate star forma - 
tion in general dGoodwin et al.ll2007t iDuchene & Krausll2013l) . 
Because binary properties reflect the main characteristics of bi- 
nary formation, they may help us determining the most com- 
mon mechanisms for the formation of multiple stars. In a bi- 
nary system of stars with masses Mi and M2 (M\ > M2), the 
mass. ratio is conventionally defined as q — Mi/M\, Similar to 
the initial mass function (IMF) for single objects, the compan- 
ion mass-ratio distribution (CMRD) is the distribution of q val- 
ues as a function of primary mass. Tidal capture models predict 
that for each primary star the mass of the secondary is chosen 
randomly from the single-star mass functio n, and the CMRD re - 
flects the IMF (McDon akT& Clarkdll993l:lKroupa et alJl2003h . 
In fragmentation scenarios subsequent continued accretion onto 
both objects from a common reservoir tends to equaliz e the 
masse s, resulting in a q distribution peaked toward unity (Bate 
2000). Capture is unlikely to be the most relevant binary for- 
mation mechanism, but it may still occur during the dissolu- 
tion phase of star clusters, causing difference s in the shape of 
the CMRD as a function of or bital separation (Moeckel & Bate 
l2010tlMoeckel & Clarkdl201 ll) . Motivated by the fact that every 
theoretical model predicts a different shape of the mass-ratio dis- 
tribution and of dependency of the CMRD on the primary mass, 
we used Monte-Carlo simulations to compare the CMRD for dif- 
ferent sampl es and to study the relation ship between the IMF and 
the CMRD ( iReggiani & Mevej d2011l) . hereafter, RM11). This 
research note represents a follow-up to RM11, in which we re- 
analyze the "universal" CMRD by adopting a different statistical 
approach (SectionFJ) and some new results on the CMRD on the 
basis of recent datasets (Section|3). 



2. Universal companion mass-ratio distribution 

The CMRD appears to be uni versal over a wide range of q 
values and primary masses (e.g. lMetchev & Hillenbrand 2009). 
According to RM11, the CMRD follows a single-slope power 
law dN/dq oc qP over the separation range 1 -2400 AU and pri- 
mary mass range 0.25-6.5 M , and there is no evidence for vari- 
ation of the CMRD with orbital separation. 

In previous work w e combined samples of M dwarfs 
dFischer & Marcvl \\99% and G stars dMetchev & Hillenbrand! 
2009 ) in the field and in termediate mass stars in ScoOB2 
( Kouwenh oven et aD 120051) adopting a x 1 fit of the combined 
binned distribution to derive the power law slope, obtaining /3 = 
-0.50 ±0.29 (Regg iani & Meverl201 Ih . The choice of the statis- 
tical method was motivated by the need of c omparing our results 
with previous studies of the CMRD (e.g. Kouwenhove n et alj 
120051: iMetchev & Hillenbrand! l2009h . However, the x 2 fit of 
a binned distribution can lead to a biased estimate, in par- 
ticular for small samples. A more robust analysis is instead 
achieved through a ma ximum-likelihood estimation method 
dFeigelson & Babul 1201 lh . This approach gives a new best-fit 
power law dN/dq oc qP, with ft - -0.18 + 0.33 to the data de- 
scribed in RM11. Although the two values are consistent with 
each other within the errors, the new estimate is flatter than pre- 
viously thought. 



3. Updates to the CMRD in the field 

Rece ntly, two new s tudies of the CMRD for solar-type 
(Raghava n et al.l 1201 Oh and M-dwarf primaries (Janso net al.l 
2012) in the field have been carried out. Since they represent the 
most complete samples to date for sun-like stars and M dwarfs, 
respectively, we applied the same statistical analysis as was pre- 
sented in RM1 1 to follow up this preliminary work. 

In the first study (Raghavan et al. 2010), roughly 200 bina- 
ries with primary masses between 0.5-3 M were considered to 
determine the CMRD over a wide range of separations (10 -1 - 
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10 5 AU) and mass ratios (0.02-1). The new CMRD appears to 
be more peaked toward unity than previously observed and the 
period distribution is unimodal and roughly log-normal with a 
peak at around 50 AU. Following the methodology described 
in RM11, we used a KS test to compare the newly obs erved 
CMRD with the CMRD bv lMetchev & Hillenbran d (2009), over 
the common range of mass ratios (0.02-1) and separations (28- 
1590 AU). The KS test returns a probability of ~ 30%, there- 
fore we cannot reject the hypothesis that the data were drawn 
from the same parent population. However, when we compare 
the two samples over the common range of mass ratios, ir- 
respective of separation, the probability is only ~ 1%, point- 
ing toward a change of the CMRD with orbital radius, because 
that of Raghavan et al. (2010) covers a wider range than that of 
iMetchev & Hillenbrand (2009). We therefore tested the possibil- 
ity of a variation of the CMRD with angular separation. To do 
this we considered break points in the angular separation distri- 
bution between 10" 1 and 10 5 AU and used a KS test for each 
of them to determine the probability that the CMRD inside the 
break point is consistent with the CMRD outside. Because we 
found probabilities greater than 1 % for any possible choice of 
break point, we conclude that we have no strong evidence for a 
dependence of the CMRD on angular separation. Moreover, be- 
cause we do not expect to see r andom pairing from cluste r dis- 
solution models inside 10 4 AU dKouwenhoven et aD l2010) and 
these widest binaries are relatively rare, we need larger samples 
in the future to test thes e models. 

The second study ( Jans on et al.l 120121) consists of 85 sys- 
tems with primary masses between 0.15-0.5 M , separations in 
the range 3-227 AU and mass ratios between 0.1 and 1. For 
M dwarfs, the CMRD appears to be fiat and the period dis- 
tribution is narrower and peaks at lower values than for solar 
type primaries, indi cating a continu ous transition from higher- 
to lower-mass stars (Burgass er et al.ll2007b . In this case as well, 
we tested the newly ob served CMRD with the CMRD from 
iFischer & Marcvl ([T992) over the common range of mass ratios 
and separations. With a probability of ~56% the KS test does not 
allow us to reject the hypothesis that the newer data were drawn 
from the same parent sample. Finally, we used the same proce- 
dure as we adopted for sun-like stars, but in the range 3-227 AU, 
to explore the dependence of the observed CMRD on separation. 
We saw no evidence of this dependence either for this sample. 

Moreover, we compared the CMRD for solar-type primaries 
from Raghavan et al. (2010) w ith the new sample of M-dwarf 
primary binaries dJanson et al.l 12012ft . The KS test returned a 
probability of 30% that the two distributions are consistent with 
each other (Figure [1}. Motivated by this result and because the 
CMRD is independent of angular separation, we combined the 
two CMRDs over the common range of mass ratios. We again 
used a maximum-likelihood method to fit the distribution and 
found a power law dN/dq cc cf , with /?=0.25±0.29. While this 
slope ft is formally consistent with the one derived in Section 
[2] (within the errors), the change in sign is significant. It is also 
worth mentioning that this fit is consistent with the mass-ratio 
distribution with power-law exponent /j=-0. 10+0.58 presented 
in a r ecent study of O-type spectroscopic binaries (ISana et alj 
2012), whereas the observed CMRD for brown dwarfs (fi ~1.5) 
points toward a different formation mechanisms for these objects 
ilGoodwin 2013)). 



4. Summary 

In this research note we have presented some updates to the 
study of RM11. First, we adopted a maximum-likelihood es- 
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Sample of G stars from the field (Raghavan et al. 2010) 

/ /^/X// Sample of M dwarfs from the field (Janson et al. 2012) 
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Fig. 1. Comparison between the observed CMRDs for solar- type pri- 
maries and M-dwarf prim aries in the fie l d. The open histogram rep- 
resents the CMRD from Raghavan et al. (2010), whereas the hashed 
histogram represents the distribution from Janso n et al.l(l2012h . The KS 
test returns a probability of 30% that the two distributions are drawn 
from the same parent sample. 



timation method to re-derive the field CM RD power law, 
based on the combination of s amples (Fischer & M arcvl 19921: 
IMetchev & Hillenbrand! l2009t iKouwenhoven et alJ l2005h de- 
scribed in RM11, to show how the dependence on the bin size 
can bias the result. 

Secondly, we analyzed recent binarity studies from the field 
(Raghava n et al.l 12010: 1 Janson et al] 1201 2|) ad opting the same 
methodology as in RM1 1. The new results from Raghavan et al. 
(2010) appear to be consistent with Mete hev & Hillenbrand 
(2009) over the common range of mass ratios and angu- 
lar separations. Th e recent updates on the M-dwarf CMRD 
dJanson et al.ll2012l) are also consistent with past results. The KS 
test does not allow us to reject the hypothesis that the CMRDs 
from Raghavan et al. (2010) and Janson et al. (2012) are drawn 
from the same parent sample. In both studies we uncovered 
no evidence for a dependence of the CMRD on separation. 
Therefore we combined the two distributions and obtained a new 
maximum-likelihood fit to the field CMRD dN/dq <x qP, with 
/?=0.25+0.29. 

Since the CMRD appears to b e independent of separat ion 
and dynamical evolution (see also Parker & Reggiani 2013), it 
represents a measurable parameter of binary stars to focus on 
when investigating binary formation mechanisms. However, we 
need larger samples to look for subtle variations of the CMRD 
with separation. In the future we aim to study the CMRD in other 
star-forming regions (e.g. the ONC) and test its dependence on 
separation for wide systems. 
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